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Abstract. Interactions between a quantum system and its environment at low 
temperatures can lead to apparent violations of thermal laws for the system. The 
source of these violations is the coupling between system and environment, which 
prevents the system from entering into a thermal state, a state of Gibb's form. On 
the other hand, for two-state systems, we show that one can define an effective 
temperature, placing the system into a 'pseudo-thermal' state where effective 
thermal laws can be applied. We then numerically explore these assertions for an 
n-state environment inspired by the spin-boson environment. 
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Introduction - - Thermodynamics of small quantum objects, such as two- 
level systems or harmonic oscillators, coupled to large environments in the low 
temperature limit has recently attracted much attention [TJ [21 131 S] m response to 
the unprecedented experimental manipulation of such systems. These studies of the 
exchange between work and heat in quantum systems is commonly referred to as 
quantum thermodynamics [5] [6] [7] and asks what happens to work and heat as 
temperature nears zero; are they still well defined, and do the laws of thermodynamics 
hold? Many effects arise due to the coupling energy not being the smallest energy scale, 
as is the case in classical thermodynamics. As a matter of fact, the strong interaction 
(where we use the term strong to indicate that the weak interaction limit no longer 
applies) between the system and the bath prevents the system from entering into a 
thermal state (a state of Gibb's form), which for T = implies entanglement but 
for finite T is a more subtle issue [5J [5J QUI El HH US]- Therefore, no well-defined 
local temperature exists associated with the system itself, and naively applying the 
equilibrium temperature of the total system may result in apparent violations of the 
second law of thermodynamics [H] [15] in the form of the Clausius Inequality (CI) 
[16] 0], and violations of Landauer's Principle (LP) \17\ [3]. While there has been 
previous work on this topic ]18[ \19\ 120] [2"T] [22] which provides a solution using a 
Hclmholtz free energy calculated from the total free energy minus the free energy of 
the bath in absence of the system, we specifically address a system commonly used 
as a model in low energy physics and take a different view point from previous work 
keeping the Helmholtz free energy calculated for the system alone, but in an effective 
framework. We investigate the quantum thermodynamics of a two-level system (or 
spin-1/2) coupled to an environment, which has many applications from dissipative 
quantum mechanics to limitations in quantum computing [3] 1231 119] , such as heat 
dissipation during quantum algorithms. 

We show that at any temperature the spin reduced density matrix can be thought 
of as an uncoupled two- level system at an effective temperature T* . Using this effective 
framework, we check that the second law of thermodynamics involving T* will hold, 
the CI becomes in fact a strict equality and LP is satisfied. The idea of an effective local 
temperature was also mentioned in |24j where they redefine the quantum quantities 
of work and heat in terms of whether energy transfers change the entropy of the 
system, while our effective framework utilizes definitions that have been carried over 
from classical thermodynamics. We also show that as temperature increases, our bath 
behaves as a purely thermal bath (with T* — > T as the thermal energy scale nears 
the coupling energy), where the von Neumann entropy essentially converges to the 
thermal entropy. 

Model — The system under investigation is a two-level qubit coupled to a bath 
in the quantum limit (close to absolute zero) modeled by a Hamiltonian of the form 
H = Hs + Hb + Hi n t where it has been divided into system, bath, and interaction 
components. We start by assuming that the total system (composed of the two-level 
system described by a Hamiltonian of Hs = h(A&x +^z) and the thermal reservoir) 
is in a thermal state at temperature T =1/(3 (we set h = ks = 1 and for simplicity), 

ptot oc exp(-/3iT). (1) 

All information about the system itself is contained in the reduced density matrix 
with the bath degrees of freedom traced out, p s = Ti blPtot]- Using this reduced 
density matrix the energy fluctuations can be calculated at a given temperature 
((5Hg) = (Hg) — (Hs) 2 ) and found to increase with the coupling between system 
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and environment in agreement with the work done by Nagaev and Biittiker [H] . For a 
finite temperature, or zero temperature, and a given bath and interaction this reduced 
density matrix can be difficult to calculate, so we calculate with a general two-level 
system density matrix, 



l/l + Z X-iY 
2 \X + iY 1 - Z 



(2) 



Any two-level state can be represented in this way with X (Y ' ,Z) being the expectation 
values of the Pauli matrices Tr[axPs] (Tr[oYp s ], Tv[azp s ])- For simplicity, we assume 
the system is not coupled via the ay operator, allowing Y to be set to zero [25] . 

First we diagonalize the system density matrix. The eigenvalues of the reduced 
density matrix are 

X± = kl^P), P=VX* + Z\ (3) 



2 

with respective normalized eigenvectors of 



2P(P±Z) \-ZtP 



X 



(4) 



Physically we can interpret the A± as the probabilities of being in the \<fi+) or 
state, respectively. Now, we define an equivalent framework in terms of an uncoupled 
two-level system at an effective temperature T* , such that p s oc cxp(— f3*Hg) and Hg 
of the uncoupled two-level system is defined below. With this effective temperature 
we treat the system as if it were a thermodynamic entity obeying effective thermal 
laws. 

We define an effective "uncoupled" Hamiltonian in order to determine parameters 
such as the effective heat. In order to do this we assume the effective Hamiltonian is 
that of an uncoupled two level system, 

H* s = l -(A*a x +e*a z ) 1 (5) 

where A* and e* which are the effective tunneling energy and level asymmetry, 
respectively. The associated energy eigenvectors of Hg are 



2?7(77Te*) \VTe* 



(6) 



with rj = \J e* 2 + A* 2 , which correspond to effective energy eigenvalues of 

E* ± = ±\r,. (7) 

We assume these eigenvectors are also eigenvectors of a pseudo-thermal state in 
Boltzmann form allowing a comparison between these states and energies to the 
coupled case. From that comparison we find the effective parameters (A*, e*, and T*) 
to treat the system as an uncoupled two-level system in pseudo-thermal equilibrium. 

We have three unknowns to determine (A*, e*, and T*). Identifying the 
eigenvalues of the reduced density matrix (3) as weights from a Boltzmann distribution 
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for the pseudo-thermal state, the effective temperature can be written in terms of the 
other two variables, 

T* = Ve* 2 +A* 2 /m[A+/A_], (g) 

which must be greater than or equal to zero from the definitions of energy (7). Due to 
the freedom given to us by introducing the three effective framework parameters (A* , 
e*, and T*) we are allowed to make a choice regarding the expectation value of H s . 
Continuity can be promoted by choosing that the expectation value of the two-level 
system energy be the same in the normal and effective framework ((Hg) = (Hg)), 
obtaining the relation 

eZ + AX = e*Z + A*X. (9) 

It should be noted here that one can calculate the additional moments of the system 
Hamiltonian following a method similar to Jordan and Biittiker [9] and find that this 
choice does not guarantee the equality of the other moments (in general (Hg) ^ (H s n ) 
where n > 2) highlighting a weakness of this formalism previously undiscussed in the 
literature. Specifically for this case (H s ) = (e 2 + A 2 )/4 and (Hg 2 ) = (e* 2 + A* 2 )/4, 
which are not equivalent unless the system is in a completely mixed state with 
X = Z = 0. 

Another relation for the effective parameters arises from the eigenvectors. We 
require the eigenvectors of the effective uncoupled two-level system (6) to be equivalent 
to the eigenvectors of the coupled system (4) corresponding to being in the excited 
(ground) state of Hg. We accomplish this by defining the effective parameters such 
that the eigenvectors have the same direction in Hilbert space. In order to determine 
these definitions we set the ratio of the elements of \4>+) ( \4>-) ) equal to the ratio of 
elements from \(j)* + ) ( |</>L) ) and find 

A* 

r] = — tt{P — Z) — e* From </>_ and <p*_, and 
X 

(10) 

r]= -^r(P + Z) + e* From and <j>%. 

Setting these equations equal we can solve for e*, 

e* = A*|. (11) 

We can now use the three Eqs. (8,9,11) to solve for e*, A*, and T*. 

eZ + AX ^ ^_ eZ + AX v 

e x 2 + z 2 z ' A "IW ' 

Z + AXI eZ + AX ( 12 ) 



- pi n[ i±|] - Pln[i±£]' 

requiring T* > from the definition of T* (8) and knowing that eZ + AX < [26]. 
The effective temperature as a function of the physical temperature can be seen in 
Fig. 2, the effective epsilon and delta do not change much within our temperature 
scale under the assumption e< A used in the numerical analysis. Such an effective 
temperature has been defined for other quantum objects such as an oscillator [J. 
Now that we have our effective parameters we have the effective Hamiltonian of the 
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system, Hg = ^(A*&x + e*&z): with the system in the "equilibrium" described 
by ps oc exp(— /3*Hg). This effectively uncoupled two level system can be treated 
as if at a temperature T* and recover the reduced density matrix of the coupled 
system, as well as the expectation value of the energy. It should be noted that as the 
temperature T increases, the interaction between system and bath should play less 
of a role and the system's reduced density matrix will tend towards one of a thermal 
state at temperature T. This will result in the effective temperature converging with 
the physical temperature as the thermal energy scale nears the coupling energy. 

We now have an effective temperature and can examine the CI which states 
SQ < TdS, with SQ being the heat transferred, T is the temperature, and dS is the 
change in the entropy of the system. We examine the system in terms of the effective 
forms of the heat, temperature, and entropy. With the effective temperature calculated 
we can now look at the effective thermodynamic entropy We define this entropy as 
the negative derivative of the Hclmholtz free energy with respect to temperature, 
S(T*) = -dF/dT*, and the Hclmholtz free energy as F = ~T*\ogZ, with Z being 
the partition function of the system alone, Z = e~P E - + e _/3 E + , where E*_ and 
represent the ground and excited energies of the two- level system, Hg. Since a thermal 
state is diagonal this effective thermodynamic entropy is equivalent to the system von 
Neumann entropy, S = — \i In Aj [57] 



s = -±[-rMe-^- + e -^: 

= A, InAj = S u , 



(13) 



where \± = e 13 E± /[e 13 E - + e 13 E +}. Inserting the eigenvalues (3), this becomes 

(l + P) 1+p (l-P) 1 - p l. (14) 



S = In2- -In 
2 



Next we calculate the effective heat of the system. We go back to the definition of the 
heat as 6Q* = Tr[H*6p], which leads to 

SQ* = U e *5Z + A*6Xy (15) 

Now we know all the terms for the CI we can look at how they change and examine 
the validity of this effective CI. Combining the effective temperature with the change 
in entropy we find 

T*dS =^p(eZ+ AX)6P. (16) 

We can now write out the effective CI in terms of the differentials 

5Q* < T*dS, 

1 1 (17) 

-(e*SZ + A*SX) < —(eZ + AX)SP. 

Simplifying this equation and substituting back in for P, e* , and A* we end up with 

eZ + AX ^^ eZ + AX vxv 

x*Tz* z5z+ x*Tz* x5x 

(18) 

<(eZ + AX) , 
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finding that the effective CI is identically an equality, 



SQ* = T*dS. 



(19) 



When we consider the effective temperature as the correct thermodynamical quantity 
to be used with thermodynamical laws, then there are no violations. 

Numerical Simulations — To further our exploration of the thermodynamics of 
these systems we conducted numerical simulations. We modeled a qubit linearly 
coupled to a bath of harmonic oscillators with M finite modes, each with only the 
N lowest energy levels. We define our Hamiltonian according to the spin-boson model 



and impose a finite number of modes, truncated after the TV lowest levels. The 
harmonic oscillator bath is described by u>i, the frequency for the i th mode, and 
the creation and annihilation operators a- and di. The interaction is linear and 
controlled by the coupling parameters, A^. We are dealing with low temperatures 
so most of the population of the bath is contained in the lower levels giving merit 
to this approximation. We set the combined bath + system to a thermal state at 
temperature T and trace out the bath environment to get the reduced density matrix 
of the system. This allows the computation of the system expectation values. 

For our examination we set our parameters in the regime of e <C A (specifically 
e = 0.01 A for the results presented), and varied the temperature between T = and 
T = A as well as the coupling constant, a. We define a such that Aj = \/2au>i, 
this definition of a is inspired by the spectral density function of the Spin Boson 
model 28 but modified to account for a bath of finite discrete modes in such a way 
that when we take the number of modes to infinity we recover the Spin Boson density. 
We investigate the simple case of a single mode bath and truncate to only two energy 
levels in the mode (more complex baths exhibit similar behavior as this basic case). 
All of the information needed to examine the system is contained in the X, Y, and 
Z expectation values (for our system Y = [25]). We find that as the temperature 
increases the system becomes more mixed and the expectation values tend towards 
zero, representative of becoming a fully mixed state as expected for thermal energies 
approximately equal to the system energy scales, as expected from previous work [29] . 

With these expectation values we can examine the entropy of the system, both the 
von Neumann entropy as well as the thermodynamic entropy based on the Hclmholtz 
free energy as before, but using the actual temperature and energy as opposed to the 
effective ones. The von Neumann entropy is defined as in Eq (14). Plots of these two 
entropies can be seen in Fig. 1. 

Mapping the spin boson model to the anisotropic Kondo model [29] , we find that 
the entropies coincide when T m Tk (the Kondo temperature) and the reservoir acts 
as a thermal bath. This is a similar result as found by Horhammer and Biittner [3] 
for a harmonic oscillator chain. The deviation between the entropies near T = is 
a result of the interaction between system and bath, which entangles the system and 
bath causing the system to be in a mixed state, and therefore S v ^ 0. While S(T) — > 
at T — 0, which can be understood by analogy to the anti-ferromagnetic anisotropic 
Kondo model [29l [28j [30] where the spin is "fully screened" at low temperatures. This 
leads to violations of LP as can be seen in Fig. 2. 



[Ml 1301 US] 




(20) 
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Figure 1. (Color Online) Thermal entropy, S(T) and von Neumann entropy, S„ 
as a function of temperature, T, in units of A. The purple horizontal line signifies 
the maximum entropy of a two-level system (S = log 2). For this data a = 0.1. 




Figure 2. (Color Online) The differences between the LP seen normally and in 
the effective framework. One can see that there is an apparent violation when 
looking at the change in heat over change in entropy, 6Q/dS, but switching to the 
effective heat results in no violation. Note that T* converges to T as T increases. 



Landauer Principle — As in the work of Horhammcr and Biittner [3] we can 
look at how this quantum two- level system relates to the LP [T7]. We can define our 
differential of heat as SQ = TdS(T) and look at the ratio 8Q/dS ul where we vary 
these parameters with respect to temperature. The LP states that the ratio of heat 
released to change in entropy should be greater than or equal to T, 

5Q/dS v > T. (21) 

Physically the increase in entropy is interpreted as erasure of information. In our 
system, full erasure corresponds to a change in entropy of log 2 (from Eq. (14)). We 
can plot this ratio for our two-level system coupled to a thermal bath and compare it 
to this inequality (Fig. 2). From this plot we can see that at low temperatures there 
is a violation of the LP, very similar in behavior to what was seen for the harmonic 
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Figure 3. (Color Online) The differences between the LP seen normally and in 
the effective framework. This shows the violation as a function of the parameter 
a. The temperature is set at T = 0.05A for this plot. 



oscillator system [3] showing that this behavior may be a general result for interacting 
quantum systems. 

Following the procedure put forth by Kim and Mahler [3] and discussed above we 
can define an effective heat and look at the LP in terms of these new quantities (Fig. 
2). The LP ends up being an equality in this effective framework, as suggested by 
the analytical solution. This approach can restore the validity of the LP, but only at 
the cost of introducing an elevated effective temperature, originating from the strong 
environmental interaction. 

Conclusion — We have shown that for a two-level system coupled to any bath one 
can define an effective temperature and place the system into a pseudo-thermal state. 
At low temperatures, the strong interaction between the system and environment 
prevents the system from entering into a thermal state leading to violations of thermal 
laws for the system, though the overall state is thermal (even at T = 0). Even 
though there are some weaknesses to this effective framework (higher moments of the 
Hamiltonian not being equivalent), once we define this pseudo-thermal state there are 
no violations of the CI or LP and the system should follow according to the effective 
laws of thermodynamics. At high temperatures |29j . this effective framework converges 
with the classical thermodynamics. 
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